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In this work we study several models of decoherence and how different quantum maps and algorithms react 
when perturbed by them. Following closely Ref. 1 1], generalizations of the three paradigmatic one single qubit 
quantum channels (these are the depolarizing channel, the phase damping channel and the amplitude damping 
channel) for the case of an arbitrarily-sized finite-dimensional Hilbert space are presented, as well as other types 
of noise in phase space. More specifically, Grover's search algorithm's response to decoherence is analyzed; 
together with those of a family of quantum versions of chaotic and regular classical maps (the baker's map 
and the cat maps). A relationship between how sensitive to decoherence a quantum map is and the degree of 
complexity in the dynamics of its associated classical counterpart is observed; resulting in a clear tendency to 
react the more decoherently the more complex the associated classical dynamics is. 

PACS numbers: 03.65.Ca, 03.65.Yz, 03.67.Lx 



I. INTRODUCTION 

The fundamental dualism inherent to quantum mechanics' 
axiomatic base, (deterministic) unitary evolution and (proba- 
bilistic) measurement, has been disturbing physicists (among 
others) ever since the very conception of quantum mechanics 
itself 1 2]. This very fundamental problem is in turn related 
to a potential practical application which has become a major 
subject of research in the last two decades: quantum infor- 
mation processing, for which the maintenance of the relevant 
system's coherence is a necessary condition. Therefore, un- 
derstanding the quantum-to-classical transition and, in partic- 
ular, decoherence (loss of quantum information) has become 
of central interest |3|. 

Decoherence's basic idea 0Sl consists of assuming that a 
closed system's evolution is always unitary and that the sys- 
tem of interest (e.g. a quantum computer), which from now on 
will be called the principal system, is never isolated from the 
environment. In this scheme the composite system formed by 
the principal system and its surrounding environment (namely, 
the rest of the universe) is a closed system and thus evolves 
unitarily. But if the accessible information is only that corre- 
sponding to the principal system, or if the access to the en- 
vironment's information is just an averaged one, by tracing 
out the environment's degrees of freedom one obtains for the 
reduced density operator of the principal system, p, the fol- 
lowing evolution 

p' = ${p) = J2^pMl , (1) 

where p' is the transformed reduced density operator |6] of 
the principal system after a certain time t. The right hand 
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side of equation Q is the so called Kraus representation (or 
operator sum representation) of a superoperator $, where M M 
are arbitrary operators satisfying the trace preserving condi- 
tion MjtM^ = I- Equation Q also guarantees that $ be a 
completely positive map and that the evolution be Markovian 
0] . A superoperator is thus a completely positive trace pre- 
serving linear map (CPTPLM), and the most general marko- 
vian evolution is ruled by such a map (the unitary evolution is 
just a particular case, when there is only one Kraus operator 
in the sum). During the rest of the paper we will stick to the 
following notation for $: 

N^M,0Mt , (2) 

where the symbol symbol "0" is defined by comparison with 

In Ref. llj] generalizations of the depolarizing and phase 
damping channels to the case of an arbitrarily-sized finite- 
dimensional Hilbert space were introduced. They were stud- 
ied in the context of the chord representation 1 8] and in a (dis- 
crete) phase-space-based approach. In particular it was found 
that, using the symplectic invariance of the chord operators, 
special decoherence models could be constructed that pro- 
duced decoherence towards selected "pointer states" by dif- 
fusing on phase space lines. 

Here we will present a brief revision of these noise models 
and include the generalization to the case of amplitude damp- 
ing channels. The focus of the present paper, in contrast to 
what was done in 1 1 ], will be to study the effect of these noise 
models on the evolution of otherwise unitary maps. We study 
the spectra of the superoperators and use DPS representations 
to display the noisy evolution. 

More specifically, we study various quantized maps that 
have different classical limits with the purpose of understand- 
ing how their underlying regularity or hyperbolicity react to 
the noise. We also study the Grover's search algorithm, which 
has no classical analog but is one of the simplest examples 
of quantum algorithms that can outperform their classical 
analogs. 
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The paper is structured as follows. In Sec. Ill Al and III B I we 
review the generalized depolarizing and phase damping chan- 
nels, show their actions in the DPS and introduce quantum 
circuits to implement them. Sec. IHCl is devoted to the gener- 
alization of the amplitude damping channel. After that, in Sec. 
Ill Dl a reinterpretation of these channels' actions is provided 
in terms of their spectra. The response of unitary maps when 
decoherence is introduced is studied in Sec. [HI] In particular, 
in Sec. IIII Al we study the action in DPS and the spectra of 
these noisy channels composed with the unitary Grover trans- 
formation; and in Sec. 1111 Bltltc composition is made with the 
quantum baker's map and with the quantum cat maps. For 
these cases, the loss of purity of quantum states is quantified 
by means of the linear entropy. Finally, we conclude with our 
results' implications in Sec. IIVI 
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H. GENERALIZED NOISY CHANNELS 



Quantum information processing's non trivial advantages 
over the classical arise when considering a large Hilbert space. 
Thus in order to study the noisy evolution of a non triv- 
ial quantum information processor, first thing we need to do 
is develop generalizations of the one-single-qubit quantum 
channels already mentioned to the case of an arbitrarily-sized 
finite-dimensional Hilbert space. 

A. Generalized depolarizing channel. 

In (l| such a generalization was proposed for the depolar- 
izing channel. The Kraus representation of the superoperator 
obtained therein in the chord representation is: 
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Figure 1 : Wigner function representation of the evolution of an initial su- 
perposition of two coherent states under the action of $dc- The vertical 
and horizontal axes correspond to the coordinates p and q respectively, and 
the gray scale intensity is proportional to the value of the Wigner function in 
that point, being white the most negative value and black the most positive 
one. In the first graphic both the interference fringes due to the periodicity 
of the DPS and the quantum interference fringes (the ones present only in 
coherent superpositions) can be seen; after applications of $£>C these inter- 
ferences disappear and, in the final graphic, they are completely replaced by 
the chessboard-like Wigner function corresponding to the completely mixed 
state jt , the microcanonical distribution. 



= (l-e)iQi + e$ DC , (3) 

where the subindex a is a shorthand notation for one of the 
TV 2 DPS points (q,p). The map $dc is also a properly nor- 
malized superoperator and its physical interpretation is quite 
simple from a phase space point of view. It performs all pos- 
sible displacements with equal weight thus averaging over all 
DPS points. It can therefore be considered within the family 
of the "diffusive superoperators" II 2u1 3Ul4l but with a very 
particular type of diffusion, one in which the density matrix 
is spread over the whole DPS uniformly. So, with probability 
1 — e, $ oc leaves the state unchanged while, with probability 
e, it averages it over all DPS points. 

In Fig. \l\the action of $dc on a quantum state can be 
seen in phase space. There we have plotted the discrete 
Wigner function of the initial "cat state" \tp) — -^(\(q,p) — 

(0.4, 0.25))+|(g,p) = (0.6, 0.75))), where | (q,p)) represents 
the coherent state centered at the phase space point (q, p), after 
no application of the map, po, and after one, two and three ap- 
plications, pi, P2 and p3, respectively. The size of the Hilbert 



space is N — 32 (five qubits) and the decoherence param- 
eter used is e = 0.8. In the first plot we can observe the 
two gaussian black spots at the centers of the two coherent 
states, the quantum interferences between the two coherent 
states (like that observed right in between the two black spots) 
and the interferences with the images that come from the pe- 
riodic boundary conditions chosen in the discretization of the 
phase space (which can now be thought of as a torus). Af- 
ter three applications of the map all interferences have disap- 
peared and the resulting graphic is the chessboard-like Wigner 
function that corresponds to the state -k, the microcanonical 
distribution, the minimal (classical) information situation. So 
not only does $£>c provoke loss quantum information (deco- 
herence) but also of the classical one. The same result is ob- 
served when applying the map to position or momentum states 
superpositions. Thus we see that the generalized depolarizing 
channel has the completely mixed state as its only invariant 
sate (as is also the case for the one-single-qubit channel). 

In view of the latter, it is now easy for us to think of a circuit 
that implements the action of this superoperator for the case 
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of N — 2™ (that is, for a set of n qubits). All we need is a 
circuit that leaves the density matrix of the n qubits unchanged 
with the correct probability (1 — e), or that replaces it by the 



state jj (with probability e) 11511 . And that is precisely what 
the circuit in Fig. [2] does. There we can observe an upper 
line marked with the symbol "/", that means transport of n 
qubits (the principal system), a central line that also transports 
n qubits (which represent the environment, initialized in the 
state jj) and a lower line that transports just one ancilla qubit 
(initialized in the state (1 - e)|0)(0| + e|l)(l|). This ancilla 
qubit works as the control for the n-qubit-controlled-SWAP 
gate, that does nothing when the state of the control is |0) 
(probability 1 — e) or exchanges the states of the two n-qubit- 
systems when the control is in state |1) (probability e). 

P -i 



(1-e) |0><0|+e|lxl| 



This superoperator is also diagonal in the chord representa- 
tion, it is: 
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(l-e)i(3l + e$ PDC . 



(5) 



Again, $pdc lS a properly normalized diffusive- 
superoperator-like map, but instead of performing a uniform 
diffusion over the whole DPS as $ p>c does, it displaces the 
density matrix to all the N points along the line of equation 
q = with the same weight -k. It diffuses the density matrix 
along the vertical direction. So $pdc does nothing with 
probability 1 — e and projects the Wigner function onto the 
vertical lines that correspond to the position states (with 
which we have arbitrarily identified the computational states) 
with probability e. But the vertical direction is not at all 
a preferential one, it only appears because of the arbitrary 
identification we have made of the computational states with 
the position states. The general case was obtained in with 
what was called the generalized phase damping channel on a 
line £rn ,7i2,»i3 



Figure 2: Circuit implementation of the generalized depolarizing channel 
for a principal system consisting of a set of n qubits (upper line). The envi- 
ronment is modeled with another set of n qubits initialized in the completely 
mixed state (central line). There is also an ancilla qubit (lower line) that works 
as the control of the n-qubit-controlled-SWAP gate. 



Generalized phase damping channel and phase damping 
channel on a line. 



The generalizations we present here were also developed 
in (jl . Written in terms of the skew projectors (or transition 
projectors) Py = \ (being |i) the i t/l -member of the com- 
putational basis, < i < N — 1 = 2" — 1), the Kraus repre- 
sentation of the generalized phase damping channel is 



t>PDC 



(l-e)/07 + e 



JV 
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Cij Pij 



QPJ, 



(4) 



where C y - must be an N x N real symmetric positive semidef- 
inite stochastic matrix. If we choose CV,- = , where 6ij is 
the Rronecker delta, we obtain a superoperator only with diag- 
onal projectors P^(from now on we will refer to generalized 
damping channel as the generalization for this choice of C y , 
other choices will be explicitly stated). In this case the action 
of the map is to do nothing to the diagonal elements of the 
density matrix it acts on and to multiply its non diagonal ele- 
ments by the factor 1 — e^exactly as the one-single-qubit phase 
damping channel does dBSl- F° r more general choices of 
Cij the action on the non diagonal elements is the same, but 
the diagonal ones are then altered. 
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where R is the number of points (q,p) in L nitn2}n3 (not nec- 
essarily N 1 9]) and the sum is performed over them. 

This superoperator's action in phase space can be appreci- 
ated in Fig. 13 where we have plotted the Wigner function 
representation of the evolution of an initial cat state for the 
case of n\ = 1, ri2 = — 1 and n-j — 0. As much as the simple 
$pcc (= $p,dc (0 „ 2 0) ) does with the position states (which 

are the eigenstates of T(o,p),P € Z), $pdc (1 _ x 0) does noth- 
ing to the state with probability 1 — e or takes it to an incoher- 
ent superposition (note how at the final plot the central quan- 
tum interference fringes have been practically removed) of its 
projections onto the eigenstates of Tj^-i), this map's pointer 
states. In a more general direction of diffusion, $pDC (n n 0) 
acts exactly as $p oc but with the computational basis being 
that of the eigenstates of TV ni „ 2 ) instead of the position basis; 
that is, it does nothing to the diagonal elements of a density 

's eigenstates and multi- 



matrix written in the basis of Ti 



(ni,n 2 ) 



ply its non diagonal ones by the factor 1 — e. Thus the pointer 
states for this model can be selected by choosing the diffusion 
line. In the most general case when the ordinate at the origin 

becomes 
plus 



12,0) 



i%3 is different from zero, the action $pdc, 
a combination of the one just described for $p dc { „ 
a unitary translation fH fl7ll . 

The circuit shown in Fig. 0] implements the generalized 
phase damping channels (if the choice for the computational 
basis is the position basis the resulting superoperator is $p p>o, 
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Figure 3: Evolution of the usual initial cat state upon application of the map 
BpDCn -ia)' ^ ne vame °f tne decoherence parameter e is the usual too, 
but now the dimension of the Hilbert space has been taken as TV = 64 (six 
qubits) for a better appreciation of this map's action. In the last panel we can 
see how the diffusion along the line p = —q has almost made the quantum 
interferences disappear, the remaining state is an incoherent superposition of 
the projections of the initial state on the eigenstates of T(—\ which are 
this superoperator's pointer states. 

if the chosen computational basis is that of the eigenstates of 
n2 ) then the resulting superoperator is %pDC (n „ ))- It 
is composed by n controlled-rotation gates R y (0) (that im- 
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Figure 4: Circuit model implementing the generalized phase damping chan- 
nel on a system of n qubits. If it is cos ® = 1 — e then the non diagonal 
elements of the principal system's transformed density matrix p' are multi- 
plied by 1 — e; while the diagonal ones remain the same. 



plement the unitary rotation R y (9) around the j/-axes in an 
angle 9 if and only if the state of the control is |1)) acting in- 
dependently on every one of the n pairs of qubits one from 
the principal system (upper lines) and one from the environ- 
ment (lower lines). The environment is initialized in the state 
|0)® n (which is a short-hand notation for the composite state 
in which all n qubits are in state |0)). The last gate is a pro- 
jective measurement in the product basis {|0) ® |0) % 

|0), |0) ® |0) ® |1), , |1) <g> |1) ® |1)}, which is 

nothing but n independent one qubit projective measurements. 
It can be easily shown fl7ll that the effect of this circuit on a 
density matrix p written in the computational basis is to leave 
the diagonal elements unchanged and to multiply the non di- 
agonal ones by the factor cos | . So, by choosing 9 such that 
cos | = 1 — e we get the desired action for the circuit imple- 
mentation of the generalized phase damping channel. 



C. Generalized amplitude damping channel 

We now turn to a generalization of the one qubit amplitude 
damping channel. This is a schematic model for the process 
of decay of a two level excited atom by means of spontaneous 
emission of a photon |4j,|5|]. Its Kraus form is: 

$adc = (I0)(0| + VT^|1)<1|) (|0)(0| + VT^|1><1|) 



+ e|O)(l|0|l)(O| . (7) 

Again, |0) and |1) are the computational states which, in this 
case, are the ground and excited states of the atom, respec- 
tively. The density matrix' elements transform according to 

( Poo poi \ _^ ( Poo + epn e )Poi \ 

V Pio Pn J V y/(l-e)p w (1 - e)p n J 

The non-diagonal elements are damped but (in contrast to 
what happened with $pdc) here the amplitude of the excited 
state is also damped and the atom eventually ends up in the 
ground state. 

We now allow the atom to have N energy levels (which will 
be used as the computational states) and the field at non-zero- 
temperature in such a way that there are photons in the field 
and the atom can absorb one of them making a transition to 
an upper energy level. In the one-single-qubit case the transi- 
tion from the excited to the ground state is carried out by an 
annihilation operator accompanied by a factor e (Kraus oper- 
ator v / e|0)(l|). Now we will need Kraus operators that take 

any computational state \i), with i = 0, 1, 2 N — 1, to the 

\i + /i) with a certain transition probability Pi : i+^. It is imme- 
diate to convince one self that the following Kraus operators 
will do: 

N-l 
t=0 
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The transition probability matrix Pi,i+^ must be real, posi- 
tive semidefinite, with all its elements lower or equal to 1 and 
greater or equal to (it must be a probability), stochastic, and 
such that Pi : i+^ = for all > i + fj, or i + fi > N — 1 
(so that there is a null probability for any state to be excited 
up beyond the \N — 1) or to decay down beyond the |0)). It 
is straightforward to show [17] that if the matrix p^j is also 
chosen to be symmetric, the resulting superoperator is self ad- 
joint. In particular we chose this matrix to be Gaussian with 

a mean width e. pij = Pij(\i — j\) oc e a? , so that the 



closer two states are in energy the greater the probability of 
transition between one another is. Thus the generalized am- 
plitude damping channel will have Kraus operators that raise 
any state a maximum of N — 1 levels (/i > 0), others that will 
lower it a minimum of N — 1 levels (/x < 0) and another that 
leaves the state unchanged (/i = 0). Its Kraus representation 
is the following, 



N-l 

t>ADC = 

H=-N+l 



(10) 



with given by l|9}. 
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Figure 5: Evolution in phase space of the usual initial state upon application 
of %adc- The central quantum interference fringes rapidly disappear and all 
other spots are stretched along the horizontal direction. 

In Fig. [5]we can see the phase space evolution of the usual 
cat state under the action of %adc for the usual parameters 
choice. We see there that already in the first iteration the 
quantum interferences have almost disappeared, showing how 
strong a decoherer this superoperator is. We can also see in 



the last panel how the surviving spots have been stretched a 
little bit along the horizontal direction. The latter can be un- 
derstood as follows, %adc picks any computational state and 
displaces it symmetrically to all its neighbours with a certain 
probability, but as we have arbitrarily chosen the computa- 
tional states to be the position states this shows in the plot as 
a spreading of the initial spots in the horizontal direction. 

For a circuit implementation of the generalized amplitude 
damping channel we need a circuit that takes the i th computa- 
tional state and sends it to the (i + n) th one with a probability 
Pij+p. A summarized scheme of such a circuit can be seen 
in Fig. [6] Once again the symbol "/" stands for transportation 
of and operation over n qubits. For example, R y l (0) repre- 
sents a controlled 2n-qubit gate that applies a different series 
of n rotations, to be described below, on the n lower qubits 
(those of the environment, all initialized in the state |0)) for 
every different state of the n upper ones (those of the princi- 
pal system) that work as the control. And the controlled-NOT 
marked with the symbol "/" represents n simple controlled- 
NOT gates each one acting independently on the n pairs of 
qubits one from the principal system and one from the envi- 
ronment. For the sake of clarity, we display in Fig. Qthe 
enclosed area of the circuit in Fig. [6] with every qubit line 
drawn explicitly. We can see there the decomposition of the 
n-qubit-pair gate Ry(0) into the n two-qubit gates Ry(6j), 
< j < n, which are not simple controlled rotations like 
those in Fig. |4] either. Rather, they are gates that apply the 
simple rotation operator around the y-axis in an angle Q- , 
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Figure 6: Circuit implementation of the generalized amplitude damping 
channel. All lines (gates) transport (operate on) n qubits. The marked block 
has been drawn in detail in Fig. Q 



R y (8j), to the environment's j th qubit when the principal 

system's j th qubit is in state \lj), where lj can take the val- 
ues or 1. Thus there are altogether 2n angles (2 for every 
one of the n two-qubit gates Ry (9 ? )) to choose in order to 
univocally specify this circuit's action. 

Writing the i th , < i < N - 1 = 2" - 1, computa- 
tional state in base 2, fa) — \h x ^2) \l n -i) ® \ln)) = 

\I1I2 ln-iln) (where the subindex 2 under the i stands for 

"binary representation of") , it is straightforward to see fPUl 
that after application of the entire circuit the state is taken to 

\h®si h®s 2 ?„_i®s n _i l n © s n ) (where "©" is the 

bit-to-bit sum that must be done adding modulo 2 every bit 
from a number with its corresponding bit from the other and 
where Sj can only take the values or 1 too) with a probability 



- , l 2 ,«2tBSlS2. 
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Figure 8: Detailed decomposition of the gate R y (8? ) enclosed in Fig.lTlin 
terms of simple controlled-rotation and NOT gates. If the control is in state 
|0), the applied rotation is in an angle 8® (R y (8j)); and when it is in state 

|1), the rotation is in an angle 01 (R y (8j)). 



Figure 7: Enclosed area of the circuit in Fig. |6| The upper (lower) lines 
are n and transport the principal system's (the environment's) qubits. The 

R y (Q? ) gate (drawn explicitly in Fig. |8} is the simple controlled-rotation 
gate Ry(8®) when it is l n = 0; and R y when it is l n = 1. 



being 



in 2 (f), 



for Sj = 



for Sj = 1 



(ID 



(12) 



So for every j there are 4 = 2 2 possible values which give 
us altogether (2 2 )™ = (2™) 2 = N 2 different combinations 
of cosines and sines in dl 11 1, that is exactly what we need to 
specify the N 2 transition probabilities pu+fj,. Thus making 
the following identification for the binary representation of the 

number p, pi = S1S2 s„_xs n , we get the desired action : 

a circuit that takes every one of the N computational states \i) 
of the principal system to a superposition of the N states \i + 
p) with probabilities pi 2 .i 2 +fj. 2 gi ven by Eq. (II It (which are in 
turn controlled choosing the 2n available rotation angles) that, 
after application of the measurement, becomes an incoherent 
superposition of them. 



D. The spectra. 



the N diagonal projectors Pa), while that of the eigenvalue 
1 — e is TV 2 — TV (corresponding to the non diagonal projec- 
tors Pij, i ^ j, which contain the coherences of the density 
matrix. This (N 2 — iV)-fold degeneracy at 1 — e depends on 
the choice = for the coefficients in ©. For random 
coefficients the degeneracy is slightly broken, as shown in the 
top right-hand plot. In the bottom right-hand graph the spec- 
trum of the generalized phase damping channel on a line for 
the choice n\ = 1, n% = and 713 = 2 is shown; there the 
degeneracy has been broken even more strongly and N out 
of the TV 2 — N that were initially at 1 — e have been spread 
along the circumference of center 1 — e and radius e in pairs 
(the degeneracy of these N eigenvalues depend on the param- 
eters ni, fi2 and 723). The bottom left-hand graph shows the 
spectrum of the generalized amplitude damping channel; here 
the previous degeneracies have been completely broken (the 
eigenvalue 1, of eigenoperator /, is simple and all the others 
are at most doubly degenerate. 

In view of the spectra we can now understand some fea- 
tures of the action of these maps in phase space. The fact that 
the microcanonical distribution is the only invariant state of 
$ dc is due to I being its only eigenoperator with eigenvalue 
1. In contrast, the invariant subspace of the generalized phase 
damping channels is TV-dimensional, spanned by the N pro- 
jectors Pa on the pointer states. On the other hand, $adc, 
when expressed in the chord basis in terms of translations, 
does not take a diagonal form as was the case for $uc and 
$dc- This is now not surprising due to its spectrum's very 
weak degeneracy 1 18[ . 



We now turn to an analysis of the spectra of these superop- 
erators, which in all cases except for the ADC can be obtained 
analytically. (See Fig. |9]) The top left-hand graph corre- 
sponds to the spectra of the generalized depolarizing channel 
and the generalized phase damping channel: in both cases the 
eigenvalues are at A = 1 and A = 1 — e, but their degeneracies 
differ. For the case of $ dc the unit eigenvalue (correspond- 
ing to the eigenoperator I) is non-degenerate and the other 
eigenvalue has degeneracy TV 2 — 1, For the case of $pdc 
the degeneracy of the eigenvalue 1 is N (corresponding to 



III. NOISY UNITARY EVOLUTION 

As was pointed out in the introduction apart from study- 
ing the decoherence models themselves the other goal of this 
work is to characterize their effects on unitary maps, to study 
the noisy unitary evolution. For this aim we will model the 
noisy evolution with two- sta ge s uperoperators as is done in 
the literature CI E E [3Tl3l- That is, given a unitary 
superoperator U, such that p' = U(p) — UpU^, we will 
consider a map composed of a unitary step plus a pure de- 
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Figure 9: Numerical spectra of the maps (clockwise) $dc an d $pdc. 
$PDC f° r a random choice of the coefficients of Cy, $l 102 ; and §ADC 
in the complex plane. The parameters are TV = 32 and e - 



= 0.4. 



coherent step described by a superoperator $ (which, for our 
case, will be one of the studied channels), $ U, such that 
p' = (% U)(p) = $(U(p)). In this approximation the total 
evolution occurs by alternating unitary and noisy steps. 



A. Noisy Grover's algorithm. 

The first unitary map we will study is Grover's search algo- 
rithm which has no classical analog but is one of the simplest 
examples of algorithms whose quantum versions are more ef- 
ficient than the classical ones |4j,|5|]. The algorithm consists 
of the successive application of the unitary operator Uq on 
the state \t/j) initialized as the uniform superposition of all N 
computational states (once again taken as the position states), 
= -A= ^^Jn 1 k) (which is equivalent to the first element 
of the momentum basis, \ip) = \p = 0)). It is designed so 
that after an optimal number of iterations T w j\fjjf the 
M marked elements \w) will have a large probability in the 
computational basis. In turn, the operator Uq can be decom- 
posed as Uq = U^Uo- Uo represents the call to an oracle 
whose only capacity is to distinguish the marked items from 
the rest. U$ is the inversion-about-the-mean operator that con- 
tains no information about the marked items and can be writ- 
ten as £/</, = I - 2\ip){ip\. 

In Fig. ^2 (top) we can see the action of the algorithm in 
phase space for N — 32, M = 1 and w = 30 (only one 
marked item, the penultimate one). Initially the state of the 



quantum computer is the momentum state \p = 0) uniformly 
distributed over all coordinate states. After every iteration 
the probability gradually concentrates on a position eigenstate 
centered on the marked \w). We can also see the growth of the 
success probability p s (initially equal to ^ = 0.0312) in such 
a way that if a measurement is performed right after the fourth 

iteration (T w j \f^i = f ~ 4) we w ^ nave more than 
a 99 per cent chance of getting \w) as the result. Note that in 
the fifth iteration this probability starts decreasing again (that 
is why it is so important to perform the measurement right at 
the optimal iteration). 

While the action of the unitary algorith in phase space was 
already studied in 13 IH; the spectrum of the map provides 
a perspective from a different point of view. It is possi- 
ble to write the initial state of the quantum computer as a 
combination of a uniform superposition of the N — M non 



marked computational states, 



J2 a eN-M \i)< 



VN-M ^q£N-M 



* 0). The 



and a uniform superposition of the M marked ones 

^E. e MM;thatis,|^ = yfZjPM- ~ 
particular thing about these two vectors is that they span a 
2-dimensional invariant subspace of Uq in which its matrix 
representation is just a 2 x 2 rotation matrix in an angle 8, 



such that sin6> 



2y / (N-M)M 



N 



This leads us to the well 



known geometrical representation of the algorithm the 
initial vector \tp) is rotated in the a/5-plane in an angle 9 per 
iteration up to the optimal iteration T when the projection of 
on |/3) (and thus the success probability p s ) is maximized. 
If the iteration process is kept on, the vector \ip) continues 
to rotate and p s starts to decrease again and so forth. The 
eigenvalues of this rotation matrix are A± = e with eigen- 
vectors | A±) = ^TjGc*) i These two eigenvalues of the 

operator Uq give rise in turn to four eigenvalues of the map 
U G (ee U G U G ): e ±2i9 and e ±i0 = 1, of eigenoperators 
|^±}(A±| and |A±)(A T | respectively, which analogously span 
a 4-dimensional invariant subspace for Uq. As far as this sub- 
space is concerned, the action of the map would be exactly the 
same regardless of the eigenvalues corresponding to operators 
outside the subspace. If the density operator is initialized in 
the state \ip){ip\ it never leaves the subspace. Therefore, the 
algorithm would work equally well no matter what the other 
eigenvalues are. We shall now see that when we consider the 
map composed by Uq and some of our channels the subspace 
spanned by |A±)(A±| and |A±)(A T | is no longer invariant. 
Thus the relationship of these four principal eigenvalues and 
all the others begins to be relevant for the search's efficiency. 

In Fig. we can see the numerical spectra of the 

maps (from left to right) Uq,($ dc )o{Ug), ($pdc)o{U g ) and 
($Adc)o(Uq) for N — 32 and e = 0.4. The first circle shows 
the eigenvalues of the unitary map by itself. The eigenvalue 
1 and those in the second and third quadrants are very degen- 
erate while those in the first and fourth ones are simple. The 
latter are the two eigenvalues corresponding to the eigenop- 
erators |A±)(A±|, while those corresponding to |A±)(A T | are 
concentrated at 1. 

In 11711 it was shown that the action of $ r>c ° n the spectrum 
of unitary maps is to contract the spectrum uniformly in the 
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Figure 10: Numerical spectra of the maps (from left to right) Uq, ($dc)°{Ug)> {$pdc)°{Ug) an d ($adc)°(Ug) m me complex plane for N = 32 and 
e = 0.4. All the eigenvalues of the unitary map Uq lie on the unitary circumference. $dc contracts the spectrum radially without breaking any degeneracy 
except for the eigenvalue 1, %poc a l so contracts it radially leaving the eigenvalues' phases intact but inducing a slight splitting of some eigenvalues. For last, 
%ADC breaks all degeneracies strongly in the radial and angular directions. See text. 



radial direction by a factor 1 — e except for the eigenvalue of 
the eigenoperator / that is left unchanged at 1. If an eigen- 
operator of the unitary map is orthogonal to / then it is also 
left intact with its associated eigenvalue reduced by the factor 
1 — e; but if it is not, then it is its orthogonal complement to / 
what is left as an eigenoperator of the composed map with the 
associated eigenvalue reduced by 1 — e too. And indeed, this 
is what can be seen in the second graph where all eigenvalues 
have been projected radially from the unit circumference to 
the circumference of radius 1 — e (also drawn explicitly) ex- 
cept for one simple eigenvalue that has been left at 1, that of 
the identity operator. As forseen before, the space spanned by 
I A±)(A± | and |A±)(A=p| is no longer invariant under the action 
the composed map, for while | A±) (A T | are traceless operators 
(and thus orthogonal to I), |A±)(A±| are not. So |A±)(A T | are 
still eigenoperators of the composed map now with eigenvalue 
1 — e, but not |A±)(A±| whose orthogonal complements to I 
are multiplied by (1 — e)e ±2te and their projections on / by 1. 
We conclude this way that in this case there is a competition 
iteration by iteration between the tendency of the algorithm 
to take the initial state to the state of interest |u>)(ui| and the 
tendency of the noise to take it to j^. 

In the third graph we can see how some degeneracies of 
those eigenvalues that were at 1 and now are around the point 
1 — e have been broken. But the most important thing is that 
the eigenvalues left at 1 are now N, those corresponding not 
to I but to any diagonal projector Pa (or alternatively, to any 
diagonal matrix). The loss of information is much less abrupt 
in this process than in the previous one, so we expect the al- 
gorithm to work better when perturbed by the phase damping 
channel than by the depolarizing channel. 

In the fourth graph we can see how %adc breaks all degen- 
eracies strongly in the radial and angular directions. This time 
there is only one eigenvalue left at 1, that of / as in the case of 
%dc- But, in contrast, there are many eigenvalues forming a 
quasi-continuum close to 1 and, what is most important, two 
of the principal eigenvalues (in the first and fourth quadrants) 
are not reduced as much as in the previous cases (they are out- 



side the (1 — e)-radius circle). So it is reasonable to expect a 
better functioning of the algorithm for this case too. 

All the conclusions drawn above from the spectra can be 
confirmed in phase space comparing the unitary (Fig. II lltop) 
and the noisy cases (Fig. ^Jbottom and Fig. I12> . As in the 
upper plot of Fig. ^2 the lower is the Wigner function rep- 
resentation of the evolution of the computer's state initialized 
as p = \ip)(ip\ = \p = 0}(p = 0| but under the action of 
the composed map {$dc)o(Ug)- We can see how all success 
probabilities are drastically reduced with respect to the uni- 
tary algorithm. But another important aspect to notice is that 
the highest success probability is now attained at the third it- 
eration. Evidently, at the fourth iteration (the former optimal 
iteration) the tendency from the noise to take the computer to 
the completely mixed state (evident from the gradual appear- 
ance of the chessboard-like pattern) becomes dominant over 
the tendency from the algorithm to find the chosen item. Fig. 
[2]top corresponds to ($pdc)o(Ug)- This time the maximal 
success probability is attained also at the third iteration instead 
of the fourth one, but the values obtained for the probabili- 
ties are about two times as high as those for the depolarizing 
channel case in accordance to what was anticipated in the dis- 
cussion about the spectra. The lower graph of Fig. H2lcorre- 
sponds to ($adc)o(Ug)- Here the probabilities are a little bit 
lower than for ($pdc)o(Ug) but certainly still much higher 
than for ($dc)o(Ug) and the optimal iteration is the fourth 
one as in the unitary case, again in accordance with what pre- 
dicted from the spectra. On the other hand the lines parallel 
to that one of the marked item observed in the last iterations 
are nothing but the effect of the horizontal spreading %adc 
provokes. 



B. Noisy cats and bakers. 

The second family of quantum maps we shall study are the 
quantum cat maps studied by Hannay and Berry 1 2 1 ] . They 
quantize the classical motion represented by the transforma- 
tions from the torus onto itself given by the symplectic matrix 
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p , p 3 =0.0312 p j, p 3 =0.25B3 p 2 , p 3 =0.6024 



p 3 , p s =0.B969 p 4 , p 3 =0.9992 p g , p s =0.B596 




p , p 3 =0.0312 p j, p 3 =0.25B3 p 2 , p 3 =0.4907 



p 3 , p 3 =0.602B p p 3 =0.t>743 p 5 , p s =0.4539 



p , p 3 =0.0312 p j, p 3 =0.1901 p a , p 3 =0.3109 



p , p 3 =0.0312 p j, p 3 =0.1941 p a , p 3 =0.366B 



p 3 , p 3 =0.327B p 4 , p 3 =0.2633 p 5 , p s =0.1702 



1 1 1 


1 1 


1 


1 


1 






Figure 11: (Above) Unitary Grover's algorithm in phase space. After only 4 
iterations the initial momentum state \p = 0} is practically transformed into 
a position state localized precisely at the position of the marked item. From 
then on the success probability p s starts decreasing again (see text). (Bottom) 
Grover's algorithm perturbed with the depolarizing channel, the parameters 
are the same as those used in Fig. 1121 



Figure 12: DPS representation of the noisy Grover's algorithm (N = 32 and 
e = 0.3). The algorithm's efficiency is altered both in a global decrease of the 
success probability p s with respect to the unitary case and in the value of the 
optimal iteration T. The algorithm is more resistant before the phase damping 
channel (top) and the amplitude damping channel (bottom) than before the 
depolarizing channel (Fig. II ll bottom). See text. 



—G^;!). <13) 

where a and (3 are integers and (q',p') are the coordinates of 
the transformed point (q,p). The dynamics implied by this 
map is determined by the eigenvalues of M. A hyperbolic 
map is obtained by choosing a = (3 = 1. In this case we ob- 
tain \h± — 2 ± y/3, with eigenvectors (q±,p±) — (1, ±\/3) 
along the stable and unstable directions. For — a = j3 = 1 we 
obtain A^ ± = ±i with no real eigenvectors, and the map for 
this case is an elliptic rotation in phase space, whose dynam- 
ics is completely regular. The parabolic case, representing a 
phase space shear, can be obtained (for example) with a = 



and /3 = 1; 

The quantum propagator for this classical maps is given by 
the unitary operator Uc, whose matrix representation in the 
discrete position basis is 

U G (q',q) = (q'\Uc\q) 

1 1 2 

= c -% L »-Fi(t/%?) _ c -% l »(q'?' -q-q+Pq 2 ) 

Vn Vn 

where N as usual stands for the dimension of the Hilbert space 
and F\{q', q) = aq> — q q + f3q 2 is the generating function 
that plays a fundamental role in the quantization procedure 
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Let us now introduce the last map we shall study, the 
baker's map. This is one of the simplest maps displaying 
strongly chaotic behaviour and, in spite of its simplicity, it 
possesses a very rich dynamics both in its classical and quan- 
tum versions. The map is is an area-preserving transformation 
defined in the [0, 1] x [0, 1] phase space square (the torus with 
periodic boundary conditions) as 

q' = 2q-[2q], 



Depolarizing . 



(15) 



where the square brackets symbolize the integer part of the 
number between them. The transformation has a very sim- 
ple geometrical interpretation, as a "stretching" step followed 
by a "cutting" step, as a baker rolling a dough. The map is 
uniformly hyperbolic with a single Lyapunov exponent = 
ln(2). Moreover, at every point the stable and unstable mani- 
folds are parallel to the coordinate axes (((?b_ , PB- ) — (lj 0) 
and {qB + ,PB + ) — (0, 1). To quantize it we follow the quan- 
tization procedure of Balasz and Voros 1 22], which yields the 
following unitary quantum propagator: 



Up = Ul 



N/2 









(16) 



where U f n is the change of basis matrix from the position to 
the momentum basis (that is, the discrete Fourier transform) 
and whose matrix elements in the position representation are 
U Flf (q>,q) = (q>\U FN \q) = ^e~^. 

We quantify the action of the noise on the map through the 
linear entropy S, defined as S = —ln(Tr(p 2 )). The minimum 
value that this quantity takes is Smui — 0, which corresponds 
to a pure state; while the maximum is Smax = ln(N) and 

corresponds to the completely mixed state -k . 

In Fig. [OJwe have plotted the evolution of the linear en- 
tropy of an initially coherent state as a function of the number 
of iterations of the three generalized channels composed with 
the four quantum maps. In all cases we have taken N = 32 
and e = 0.2. 

The uppermost plot corresponds to $_dcs the completely 
depolarizing and most degenerate superoperator. We can ob- 
serve there how it induces the same entropy growth for all four 
maps: there exists first a linear regime and then an asymp- 
totic tendency to saturation at the value Suax = ln(N) = 
ln(32) w 3.47. So, as far as S is concerned, it is as though 
the application of $p,c "erased" the dynamics imprinted on 
the state by the previous unitary step. And, in view of the dis- 
cussion of the last paragraph of Sec. Ill Al that is not much 
of a surprise; because we know this map's action can be un- 
derstood as leaving the state intact, with probability 1 — e, or 
taking it to jj, with probability e, no matter what the initial 
state was (or, in particular, no matter what the applied unitary 
map was). 

In the central plot we can see the evolution of the entropy 
but this time for the composition of $ pp,c (whose spectrum is 
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Figure 13: Linear entropy S as a function of the number of iterations for 
$DC (top), $PDC (center) and %adc (bottom) composed with the four uni- 
tary maps; N = 32 and e = 0.2. The computer has been initialized in a 
coherent state. For the case of $adc the number of iteration has been taken 
much greater so as to fully appreciate the evolution of S until it reaches sat- 
uration (note the difference in scales). The behaviour of the entropy depends 
on both the degeneracy of the decoherent superoperator and the degree of 
complexity of the associated classical map (see text). 



less degenerate than $£>c' s ) with the four maps. In this case 
the situation is slightly different, for the curves corresponding 
to {$pdc) (Ub), (§pdc)o(U h ) and {% PD c) {Up) display 
pretty much the same behaviour as before: but on the contrary, 



11 




for the case of (%pdc)o(U e) a little slower entropy growth is 
observed. 

Finally, in the lowermost graph we can see a completely 
different situation for $adc, the superoperator with the least 
degenerate spectrum we have studied. The curves correspond- 
ing to ($adc)o(Ub) and ($adc)o(Uh) are almost identical 



and, before saturation, they present a linear growth regime as 
before but with a considerably lower slope (compare the dif- 
ference in abscissa scale with the two previous cases). The 
curve corresponding to ($adc)o(Up) also possesses a lin- 
ear growth regime but with an even lower slope. And the 
curve corresponding to ($adc)o{Ue) has completely aban- 
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doned the linear regime and has the slowest entropy growth. 

In summary the depolarizing channel (whose spectrum pos- 
sesses the maximum degeneracy) is insensitive to the classi- 
cal dynamics of the map it acts on. For the phase damping 
channel ( with a slightly less degenerate spectrum), the elliptic 
cat map ( with a regular classical behaviour) shows a slightly 
lower rate of entropy production. For the amplitude damp- 
ing channel (with its spectrum barely degenerate) all four en- 
tropy production rates are much slower than the two previous 
cases; and there is also a strong dependence on the nature of 
the classical map, with chaotic maps showing faster decoher- 
ence, which could be related to the idea of associating chaos 
to a quantum system according to how sensitive it is to deco- 
herence 12311 . 

As for the case of the superoperators acting on Graver's 
transformation, here it is also possible to recover from the 
spectra the conclusions just drawn above from the entropies. 
In Fig. ^] we can see the spectra of all four unitary maps in 
the (from up to down) first row and of $dcs $p dc an d §adc 
composed with them in the second, third and fourth rows, re- 
spectively. The (from left to right) first column corresponds 
to the map Uhc, the second to U pc, the third to Uec an d the 
fourth to Ub- We can see how $dc contracts the spectra radi- 
ally leaving all eigenvalues but the 1 on the circumference of 
radius 1 — e. $ pdc does approximately the same thing with 
most eigenvalues, but, on the other hand, for all four maps 
we find some eigenvalues that are left closer to the unitary 
circumference (which means that the behaviour of the com- 
posed map is closer to the unitary one's). Finally, we can see 
a strong enhancement of this tendency for the case of %adc- 
When composed with the two hyperbolic maps we find more 
eigenvalues close to the unit circumference than in the two 
previous cases, but still concentrated at a distance 1 — e from 
the origin. When composed with the parabolic map there are 
even more eigenvalues close to the unit circumference. And 
for the elliptic map the eigenvalues with modula close to unit 
are so many that they are seen in the plot as a continuum of 
eigenvalues accumulating on the unit circumference. So, also 
from this brief analysis, we can conclude that for the elliptic 
case the composed map remains the closest to the unitary one, 
behind it the parabolic map and last the two hyperbolic ones. 

Nevertheless, we have not yet exploited all the possibili- 
ties of the decoherent tools we have developed, for in the last 
section we have not appealed to the ability of choosing the 
preferred basis of decoherence we possess through the gener- 
alized phase damping channel. For example, we can explore 
how the reaction of Ub (and other maps' too) depends on the 
direction of diffusion. And that is exactly what can be ob- 
served in Fig. ^] where we have plotted the evolution of the 
entropy of a system undergoing application of the composi- 
tion of the depolarizing, phase damping (along several differ- 
ent lines) and amplitude damping channels with the quantum 
baker's map (top), the quantum hyperbolic cat map (center) 
and the quantum parabolic cat map (bottom). 

For the case of Ub the system was initialized in a momen- 
tum state, the superoperators used were %dc> $iu , $_e,i_ 10 > 
$Loio' $£ioo an d §adc ■ We can see that, as expected, the 
curve corresponding to $ dc is located on top of all others; and 
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Figure 1 5 : Linear entropy S as a function of the number of iterations for the 
generalized depolarizing, phase damping (along several different lines) and 
amplitude damping channels with the quantum baker's map (top), the quan- 
tum hyperbolic cat map (center) and the quantum parabolic cat map (bottom). 
In all cases the parameters are N = 32 and e = 0.2. For the case of the 
baker's map the computer was initialized in the momentum state \p = 0.25), 
while for the cases of the hyperbolic and parabolic cat maps it was initialized 
in the coherent state \(q,p) = (0.25,0.25)). For the two fully chaotic cases 
there exists a slight dependence of the entropy production rate on the diffu- 
sion direction. This dependence is considerably stronger for the parabolic 
case, for which the number of iteration was taken much greater so as to fully 
appreciate the evolution of S until it reaches saturation (note the difference in 
scales). (See text). 
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that of %adc below all others. The novelty, though, is that in 
spite of all four different phase damping channels having the 
same degeneracy, there are two of them ($l 110 and $l 1 _ 10 ) 
whose curves locate right next to that of $ uc an d other two 
(Sioio an d $£100) whose curves grow a little bit more slowly. 
The only peculiarity in $l 100 and $ l ow is that the directions of 
the diffusions they perform coincide with the classical map's 
invariant directions (its stable and unstable manifolds). Equiv- 
alent results are obtained if the state is initialized in a position 
state. And also for the case of f/#c> there the initial state was 
taken to be a coherent one, and the superoperators used are 
the same except for the phase damping channels on oblique 
directions that have been replaced by phase dampings on the 
lines of equations lOOp = ±173g, whose slopes are equal 
to ±1.73 w ±V3, the slopes of the classical map's stable (+) 
and unstable (— ) manifolds. Even though these two directions 
are just approximately equal to the map's invariant directions, 
we can observe in the case of the unstable manifold a tendency 
for the entropy to remain below the others' (except for that of 
$adc, °f course). 

There seems then to be a tendency of slower entropy growth 
when the diffusion is made along an invariant direction of the 
classical map. And we can see how this tendency is consider- 
ably enhanced for the case of Upc- There the initial state is 
also a coherent one and the superoperators taken are the same 
as those for the case of Up- Note how not only does the curve 
corresponding to the line of equation p = q (the only invari- 
ant direction of the classical map) initially grows more slowly 
than that of the line p = —q (the curves corresponding to $p>c 
and to the lines of equation q = and p = have not been 
drawn for they are identical with the latter one) but also how it 
deaccelerates so much that it even becomes lower than that of 
$adc an( J finally comes to a stop at a value much lower than 
the saturation value Smax (note the difference in abscissa 
scale). So, in the same way as when we compared the maps' 
reactions before %adc (Fig- H3 bottom graph), here for the 
parabolic case one also finds a much less decoherent reaction 
than that of the two fully chaotic cases. Then, also as in that 
case, one would expect this tendency to be even more marked 
for the elliptic case; but we know that this map possesses no 
invariant directions. Moreover, for this map the entropy was 
observed to be independent of the diffusion direction and all 
curves are identical to that observed in the central plot of Fig. 
I13lfor the case of the simple phase damping channel. 

IV. CONCLUSIONS 

We have studied some superoperators that produce various 
kinds of non-unitary quantum operations that generalize the 
standard one qubit depolarizing, phase damping and ampli- 
tude damping channels to systems with arbitrary finite dimen- 
sions. The depolarizing and phase damping channels are di- 
agonal in the chord basis consisting of phase space translation 
operators and can then be thought of as a special kind of ran- 



dom unitary process. This is not the case for the amplitude 
damping channel. The efficiency by which these channels 
take a coherent superposition into different statistical mix- 
tures was shown in phase space and was also understood from 
their spectra. The generalized amplitude damping channel 
was found to possess a very weakly degenerate spectrum (in 
contrast to the completely degenerate spectrum of the depo- 
larizing channel or the also very degenerate spectrum of the 
phase damping channel). 

As a first example of the action of these superoperators 
when combined with unitary maps we chose Graver's search 
algorithm. In the unitary case the algorithm possesses an in- 
variant subspace and its operation is a simple rotation in this 
subspace. We find that our noise models do not alter substan- 
tially this picture. The invariant subspaces are now approxi- 
mate as they are weakly coupled to the rest of the spectrum. 
As a result the search efficiency is altered both in a decrease of 
the success probability and in the value of the optimal iteration 
time. But, nevertheless, the results indicate relatively good 
level of resistance to this models of noise, with success prob- 
ability values ranging fromp e w 0.25 (for %dc) to p e ~ 0.5 
(for $pdc ar, d $adc); ror a decoherence parameter e = 0.3. 

Another example of the effect of these superoperators was 
given by applying them to the quantum versions of classical 
maps possessing dynamics with different levels of complex- 
ity: the baker's map and the hyperbolic, parabolic and ellip- 
tic cat maps. This time the reactions were analyzed from the 
point of view of their spectra and of the linear entropy produc- 
tion rates induced by the composed maps. The classical prop- 
erties of the unitary maps are reflected in two ways in the com- 
posed spectrum: for hyperbolic maps, level repulsion of the 
unitary map eigenvalues prevents degeneracies and therefore 
a given level of noise acts perturbatively producing a more or 
less uniform contraction of the spectrum towards a modulus 
of order 1 — e. On the other hand integrable maps show in 
general many degenerate subspaces that respond to the noise 
differently, the degeneracy being split in a non perturbative 
way and leaving many eigenvalues close to the unit circle. In 
this case entropy is produced at a slower rate. For the phase 
damping on a line acting on hyperbolic and parabolic maps we 
find a dependence of the decoherence rates with the alignment 
between the invariant manifolds and the line of diffusion. 

In more realistic models of decoherence acting on maps of 
mixed phase space the spectra may not be as simple and de- 
generate as the ones considered in this work, but still the anal- 
ysis of the spectrum of the composed action will carry impor- 
tant information about the time behaviour of the open system. 
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